Distribution of the Absorption by Chaotic States in Quantum Dots 
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The mesoscopic fluctuations of the absorption at optical transitions from a low energy regular 
state to high energy chaotic states in an aggregate of semiconductor quantum dots is studied. We 
provide a universal dependence of the distribution of the absorption coefficient on the total number 
of dots and the ratio of the level broadening to the level spacing. The distribution remain broad 
even at large broadening, and the absorption spectrum should demonstrate a strong sensitivity to 
weak magnetic field in the region of large and weak absorption. The results can also apply to the 
absorption of Rydberg atoms in strong magnetic field at the pre-threshold ionization. 
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Semiconductor quantum dots are now a subject of in- 
tensive study as potential devices with controlled optical 
properties Because of restricted geometry, their 

spectra are discrete and are modified by variation of the 
dot size. The experimental absorption spectrum demon- 
strates a good agreement with theoretical calculations for 
optical transitions between the first few size-quantized 
levels For these low energy states, fluctuations of 

the dot size give broadening to observed transitions, but 
their levels remain still well identified. Fluctuations of 
the confinement potential play no prominent role since 
the corresponding electron wave length is of the order of 
the dot size, i.e., exceeds a length scale of these irregu- 
larities. For sufficiently high energy states whose wave 
lengths amount to an atomic scale or a characteristic 
scale of fluctuations, their levels and wavefunctions are 
determined by such 'imperfections' of the confinement 
potential. The size of the dot controls only the aver- 
age interlevel spacing. The confinement part fluctuates 
from dot to dot and even for a given dot wavefunctions of 
such high energy states changes unpredictedly from level 
to level. This allows us to take the statistical description 
for the high energy chaotic states. 

In the present paper, we study the absorption coeffi- 
cient for the transitions from a given initial state into 
these high energy chaotic states. The initial state can be 
the ground state or a low energy state as well as a local 
electronic state, e.g., at the donor inside the dot. For an 
aggregate of dots, the experimental absorption spectrum 
should show featureless background with rare high peaks 
and the weak absorption regions. Our result describes 
the statistics of these fluctuations. We have found that 
the amplitude of fluctuations in the absorption depends 
strongly on the total number of dots, and also on the 
ratio of the level broadening to the level spacing. When 
the latter is small, the fluctuation of absorption spectra 
becomes highly asymmetric, demonstrating abrupt decay 
for small ones and slow power-law decay for large fluctu- 
ations. The fluctuations remain appreciable even when 



the broadening is much larger than the level spacing and 
therefore the fluctuations should be still observable for 
an ensemble of finite number of dots. 

Similarly to the statistics of energy levels [0 and wave- 
functions in chaotic systems the distribution func- 
tion for the absorption coefficient manifests the universal 
behavior which is determined only by the generic sym- 
metry of the system ||] . In the absence of the spin-orbit 
interaction, there are two universality classes present: the 
orthogonal class for time-reversal-invariant (T-invariant) 
systems (/3 = 1), and the unitary class for time-reversal- 
breaking (T-breaking) systems (f3 — 2). The crossover 
between these two universality classes will occur with the 
magnetic flux 4> through a dot around cf> ~ h/e^fg where g 
is the dimensionless conductance |J . We investigate the 
statistics of the absorption for both universality classes. 

Analogous absorption spectrum should be observed for 
Rydberg atoms in a high magnetic field where the pre- 
ionization states are considered chaotic juJ. 

Consider the absorption for a single dot at optical tran- 
sition from a given initial state |$o) with the energy Eq 
to excited states \ip n ) with energy E n , which we assume 
fully chaotic. Under the resonance, the real part of the 
polarizability a is negligible, so a(Q) becomes 

a(Q) ^ nrY,\Pno\ 2 S r (Q-Sl n ), (1) 

n 

where £2 n = E n — Eq, and P„o = (V'n|-P| ( I > o)j the tran- 
sition dipole matrix element. The explicit form of the 
transition dipole P is not needed to proceed the follow- 
ing argument. <5p is a 5- function with a finite broadening 
r, defined by 2tt5 t {E) = T/(E 2 +T 2 /4). A basic as- 
sumption here is that the level width is fixed among these 
high energy excited states. Since such broadening is dom- 
inated mainly by electron escape into the host matrix and 
an electron can leave a dot through its whole boundary, 
fluctuations of level widths should be suppressed. 

The absorption (extinction) coefficient k in this regime 
is determined by 



1 



k(Q) = yj2ir\a(n)\. 



(2) 



Reflecting random nature of tt n and chaotic wavefunc- 
tions, a(f2) and «(fi) become highly fluctuating quanti- 
ties, either from sample to sample or by slightly changing 
Q. To examine these statistical properties, we will in- 
vestigate the following the probability distribution func- 
tions: 



V(x;{a}) = (6(x-a(Cl)/a )), 
V(x;{K}) = {5(x- K (n)/ Ko )), 



(3a) 
(3b) 



where ao and kq are the average values of a and k, re- 
spectively (see Eq.(Q) below). We obtain the analyti- 
cal expressions of V(x;{k}) and V(x;{a}), both for T- 
invariant and T-breaking systems. 

We proceed this task by connecting the distribution 
functions of a or k with that of the local density of states 
v(r) defined by 



u(r,E) = Y,\i> n (r)\ 2 5 r (E - E n ), 

n 

V{x;{v}) = {8{x-v{r,E)/v )), 



(4) 
(5) 



where — (v(r,E)}. By utilizing statistical properties 
of chaotic wavefunctions, we can show that 



V{x;{a})=V{x;{v})=V{x), 
V(x;{k}) = 2xV(x 2 ). 



(6a) 
(6b) 



First we show how Eq. (g) was derive. The average of 
the polarizability become independent of f2 and is given 
by 



«0=f (|Pn0| 2 ) 



(7) 



where A = 1/voV is the mean level spacing with the 
dot volume V. ( |-P n o| 2 ) is the average intensity of the 
transition dipole moment |llf| defined by 



(\Pno\ 2 ) = y /aVidr 2 ($o|P|ri>/(r 12 Xr 2 |P|$o>, (8) 

with /(ri 2 ) = V (^ n (T*i)V'n( T *2)) to describe spatial cor- 
relations of the wavefunction amplitude. To calculate the 
distribution, we need its higher moments 

( I^n0| 2fc ) = / I] ( dT M ( $ 0l P \Ti) (rj| P |$„; 

x(Mri)r n (r' 1 )---MrkW n (r' k )). (9) 

The spatial correlation between the local values of wave- 
function exists only within the mean free path £. Beyond 
such distance, the wavefunction can be seen to be fluctu- 
ating independently. As a result, in the leading order of 
£\ d ~ 1 /V (A is the wave length), we get 



<|^o| 2fc H(|^o| 2 ) & 



fc! 



(for 13 = 2) 



(2k- 1)!! (for (3 = 1) 



(10) 



This shows that, as well as If/ViMI 2 , the distribution of 
|P„o| 2 and equally oscillator strength f n = 2mfl n \P n o\ 2 
are characterized by the the Porter-Thomas distribu- 
tion after an appropriate rescaling. By comparing 
Eqs. (0) and (g), we can conclude Eq. ([3a]), then Eq. ( pb| ) 
follows from Eq. (|^). 

To evaluate V(x) analytically, it is convenient to work 
on its Laplace transformation C(s) flTJf 



C{s) 



dx e 



V(x), 



(11) 



where 7 = 71T/A. The averaging for C(s) can be de- 
composed into averaging over eigenfunctions and energy 
levels, and the former is recasted by the Porter-Thomas 
distribution. Since the Porter-Thomas distribution de- 
pends on the symmetry parameter 0, the result of the 
integration over eigenfunctions C(s) = Cp(s) becomes 



fdct[(E-H) 



V4]\ 



0/2 



\&et[{E-H) 2 + f 2 /A\) 



(12) 



where H is the Hamiltonian of the system and T = 
r-^/i + 4s//3. (We introduce 7 = nf/A for later use.) 

Eq. (E2) can be evaluated by the supermatrix Q 
method mM. A trick is needed for /3 = 1 to generate 
[det (• • •)] in the numerator from the Grassmann inte- 
gration. This can be achieved by using the form 



A, 00 



det[(£-iJ) 2 +r 2 /4] 



'det[(E - H) 2 + T 2 /4] \{E - H) 2 + f 2 /4] , 



which can be readily expressed within the Q-matrix for- 
malism [[TJJ. After completing the mapping, Eq. ( |l2| ) is 
found to be equal to 



JdQ exp {STr [(a A + a x A s 



a 2 C B )Q}} (13) 



where A = diag(l 4 , -I4), A s = diag(l 2 , -I2, —I2, 12), 
and Cb = diag(02, <7i, O2, —01) (<jj are 2x2 Pauli matri- 
ces) . The definitions of STr and integration of Q as well 
as the explicit structure of Q-matrix can be found in (jfJJ . 
The coefficients (i = 1, 2, 3) are for the unitary case, 

«o = (7 + 7)/8; a x = (7-7)/8; a 2 = 0, (14) 

and for the orthogonal, 

ao = (3 7 + 7 )/16; 2d = -a 2 = (7 - 7 )/8. (15) 

Evaluating C(s) and V(x) for the unitary universal- 
ity class (/? = 2) was already done in the framework of 
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the local density of states distribution |13|] . Their results 
read 



A, (s)= |r- [(7 + 7)V - (7 - 7) 2 e" 7 ] 



(16) 



7 



• exp [— 7 (x + 1 /x) /2] 



87TX 3 

x [2 cosh 7 + (x + 1/x — I/7) sinh7] . (17) 

By using Eqs. (^]), we have the distribution of the ab- 
sorption coefficient for the unitary case. 

To evaluate V Q {x) or C (s) is much more laborious but 
still durable. Technically the difficulty results from ex- 
panding the exponent of Eq. (13) and taking the highest 
order Grassmann term to complete the integration dQ. 
After straightforward but lengthy calculations, we found 
that C (s) is given by 



C ( S ) = e^-^/ 2 +2 / e -i(T+7)A 1 A 2+7 A 



[F(Ai,A 2 ,A)] 2 dAidA 2 dA 
(A 2 + A 2 + A 2 -2AAiA 2 -l) 2 



, (1* 



where the integral region are defined by Ai, A 2 € (1, 00) 
and A € (—1, 1), and with \ii — yA| — 1 (for i — 1,2), 



F(Ai,A 2 ,A) = 



7-7 
An 



dp<p(p)e^^ p ^~"< )/4 , (19) 



(pip) = (AiA 2 + - A)/ vi -p 2 . 



(20) 



V {x) is obtained by completing the inverse Laplace 
transformation of Eq. (fill). Evaluating the integral by 
deforming the integration contour leads to 



Po(x) = 



47TX 3 

where the integral kernel K is 

K(Xx, X 2 , X,pi t p 2 ) 



i ^|e-i(- 1 ) 2 + JdXjldXidpi e-^^-V-^-^/^KiX^X^X^^) 



7 . . , , 7C 3 - 27XC 2 - 6x( + 7X 2 C + Ax 2 



toPga-rwrvi ( A 2 + a 2 + A 2 - 2AA!A 2 - l) 2 
C(Ai, A 2 , X,px,p 2 ) = AiA 2 + nifi 2 (pi +p 2 )/2. 



(21) 



(22) 
(23) 



Combining Eqs. ( p7| , |2l| , p2| ) with Eq. (j^) consist of our 
main results in the paper. 

The results for T-invariant systems (the orthogonal 
class) are particularly interesting since this is the usual 
symmetry in observing absorptions in quantum dots. We 
also remark that V{x) can be observed as the local den- 
sity of states distribution in nuclear magnetic resonant 
(NMR) experiments. NMR experiments are performed 
rather within the T-invariant situation, though only the 
analytical expression was known for the unitary case so 
far jl3| . We emphasize that our obtained result Eq. j2l| ) 
serves in this respect as well. 

Next we examine various asymptotic behavior of 
V {x). For 7 ^> 1, V (x) can distribute only around 
the unity, otherwise it's suppressed exponentially. The 
dominant behavior for 1/7 < 1 < 7 is characterized by 



V (x) 



7 



Airx 3 



exp 



7 



(Vi- 



1 



(24) 



The form obtained above describes multilevel absorption. 
The number of levels which contribute to the absorption 
is of the order 7. The individual contribution are ran- 
dom because of different wave function and therefore the 
matrix elements are uncorrelated. In accordance with 
the central limit theorem, the distribution function is ob- 
tained to be Gaussian with width cx 1/^/7- 

For 7 < 1, we can decompose the behavior of V (x) 
into the three region: (1) 1 « 7 « 1/7, (2) 7 « 1 « 



I/7, and (3) 7 <C I/7 <C x. The evaluation of the asymp- 
totic behavior from the analytical expression leads to 



V (x) 




(x < 7 < I/7) 

(7 < x < I/7) (25) 

(7 < I/7 < x) 



where C is a numerical constant. Let us note that the 
distribution function Eq. (|2^) is very asymmetric. The 
maximum V (x\ {a}) — V {x; {v}) lies at x w 7 -C 1 and 
on the right side from the maximum the function decays 
by power law as x -3 / 2 in contrast with a; -5 / 2 in the uni- 
tary case. The left tail of V [x\ \y\) decay exponentially 
as exp(— 7/4x). The absorption at 7 <C 1 is determined 
by the rare single level, thereby the distribution occurs to 
be shifted to small values. The power decay is formed by 
the spectral fluctuation. At the same time far right and 
left tails is due to the fluctuation of the matrix element 
of wavefunction. 

In the final, consider the absorption from a system of 
iV-uncoupled dots which have almost identical volumes 
but different shapes. We can write down the distribution 
function of the absorption from such a system by 

/ ^ \ /•-HOC 7 

In the unitary case where C(s) is given by Eq. (|T6|), we 
can write down explicitly, 
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- — e 2 

ZTTX 



N k 
(x+l/x) * "* * 



2m- 



(cosh7) Ar_fe (sinh7) fc 



(26) 



2 k (2N~fx) 2m - k + 1 

where Hk is the Hermite polynomials and its analytical 
continuation into negative k. In Fig. 1, we present the 
distribution of the absorption coefficient vi N \x; {n}) = 
2xVu N \x 2 ) for 7 = 0.1 and 10 by changing N. When 
7^1 (Fig. la), the distribution Vi N \x) is nearly Gaus- 
sian around x = 1, but with slight asymmetry. Even for 
7 = 10 and N = 50, we see the fluctuations amounts 
to the order of 10 %. The behavior for 7 <C 1 is quite 
different (Fig. lb). For 7 <§; 1 and N — 1, there is a 
peak around x ~ 7 with strong asymmetry, and the peak 
position move gradually to 1 by increasing N. 

The A-dependence of V^ N \x) in the orthogonal case 
is qualitatively very similar to that for the unitary case, 
except for the different power-law decay of the tails. 

In conclusion, we have studied the statistical proper- 
ties of the absorption in an aggregation of A-uncouplcd 
quantum dots at the transition between a low energy reg- 
ular state and high energy chaotic states. The transition 
matrix elements and the oscillator strengths are shown 
to obey the Porter-Thomas distribution, thus the statis- 
tics of the polarizability has turned out to be identical to 
the local density states distribution. We have found that 
statistics of large deviation of the absorption from its av- 
erage value are strongly distinguished between systems 
with and without T-invariance. Therefore application of 
a weak magnetic fields should have a pronounced effect 
on the absorption spectrum by suppressing large fluctu- 
ations. 
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FIG. 1. Distribution functions Vu \ x \ { K }) of the absorp- 
tion coefficient in the unit of its average value for iV-uncoupled 
dots. Plottings are shown for (a) 7 = 0.1 and (b) 7 = 10, 
where 7 = 71T/A and V and A are the level width and the 
mean level spacing. 
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